Asymmetric field dependence of magnetoresistance in magnetic films. 
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We study an asymmetric in field magnetoresistance that is frequently observed in magnetic films 
and, in particular, the odd longitudinal voltage peaks that appear during magnetization reversal in 
ferromagnetic films, with out-of-plane magnetic anisotropy. We argue that the anomalous signals 
result from small variation of magnetization and Hall resistivity along the sample. Experimental 
data can be well described by a simple circuit model, the latter being supported by analytic and 
numerical calculations of current and electric field distribution in films with a gradual variation of 
the magnetization and Hall resistance. 
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Onsager's reciprocity relations [T] are the cornerstone 
in understanding the field symmetry of magnetotrans- 
port measurements. Magnetoresistance or longitudinal 
resistivity (measured along the current flow direction) 
is predicted to be an even function of magnetic induc- 
tion B, while transverse (Hall) resistivity is specified to 
be odd with respect to B when a magnetic field is ap- 
plied perpendicular to the sample plane. General ac- 
ceptance of these rules is so common that in numerous 
experimental cases, when current and voltage contacts 
can not be arranged in a well defined 5-probe geometry, 
the magnetoresistance and Hall effect data are respec- 
tively extracted as the even and odd in field components 
of the measured 4-probe signal. However, asymmetric 
in field magnetoresistance is quite frequently observed in 
magnetic materials (also in samples with fully symmetric 
magnetic properties and properly arranged current and 
voltage contacts) [21 [3] , although this is rarely mentioned 
and discussed [U |S] . A seeming violation of Onsager's law 
only recently attracted attention when sharp distinctive 
peaks of magnetoresistance, odd with respect to applied 
field, were found at magnetization reversal of ferromag- 
nets with an out-of plane magnetic anisotropy [51 [7] . As 
argued by Cheng et al [6^ the effect can appear when a 
domain wall (DW), located between the voltage probes, 
runs perpendicular to both magnetization and current. 
Electric fields generated by the extraordinary Hall effect 
(EHE) have opposite polarities on both sides of the DW, 
which can produce a circulating current loop and a re- 
spective extra voltage contribution. The model was used 
to explain the odd in field longitudinal voltage peaks in 
a specially designed Co-Pt multilayer film with a single 
DW gradually propagating along the sample. However, 
the effect was also observed in other samples with mul- 
tiple domains |H] , and the applicability of the " single 
wall" model in this general case is dubious. 
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FIG. 1: (a) Longitudinal voltage (o) measured in a 6 nm thick 
Ni film at 4.2 K as a function of field applied perpendicular 
to the film plane. Solid line ( ) is a fit calculated according 
to Eq. [8] Dashed lines are a guide for the eye that empha- 
sizes the high field asymmetry, (b) Hall voltage (o) measured 
simultaneously. 



In this paper we present two typical cases of asym- 
metric magnetoresistance observed in magnetic films and 
analyse their origins. We shall argue that the anoma- 
lous behaviour can consistently be explained by a grad- 
ual variation of magnetization and Hall resistivity along 
the sample. 

Fig. [l]i presents the longitudinal voltage Vi measured 
in a 6 nm thick Ni film at 4.2 K as a function of a mag- 
netic field applied perpendicular to the film plane in both 
field polarities. Thin Ni films possess the surface induced 
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out-of-plane anisotropy at low temperatures [3] respon- 
sible for the hysteresis in the magnetoresistance curve. 
Anisotropic magnetoresistance is the origin of the neg- 
ative magnetoresistance when a field is applied perpen- 
dicular to the electric current direction. Vi reaches max- 
imum at magnetic fields corresponding to the coercive 
field value when the macroscopic out of plane magnetiza- 
tion crosses zero. Notably, the magnitude of the maximal 
voltage is not equal at two field polarities although the 
location of the peaks is the same. Similar asymmetric 
maxima can be found in several publications pUlfTTlIT^ . 
The measured voltage is Ohmic (linear in electric cur- 
rent) and, if interpreted as magnetoresistance, its asym- 
metry would mean a violation of the Onsager rule. 

Fig. [T|d presents the transverse (Hall) voltage for the 
same sample. Hall voltage in magnetic films depends on 
magnetization as [I^ : 



Vt = -(RoB + iioReheM) 



t 



(1) 



where / is electrical current, t - thickness, i?o and 
Re HE are the ordinary and extraordinary Hall coeffi- 
cients, B and M are the out of plane components of mag- 
netic field induction, and magnetization respectively. A 
clear hysteresis loop is seen in Fig. [T]d which is propor- 
tional to magnetization (the EHE term contribution is 
much larger than the ordinary one; therefore we neglect 
the ordinary Hall component in the following discussion) . 

Another striking example of asymmetric in field 
magneto-voltage is presented in Fig. |2ji. The longitudi- 
nal voltage measured in a Co/Pd multilayer sample (10 
bilayers of 0.2 nm thick Co and 1.1 nm thick Pd, total 
thickness 13 nm) is shown as a function of a magnetic 
field normal to the film. The sample was prepared by 
sequential e-beam deposition of Co and Pd layers on a 
GaAs substrate. It has the six-contacts Hall bar geome- 
try, 5 mm wide and 15 mm long. The distance between 
longitudinal and transverse voltage contacts is 5 mm. 
The sample has a strong out-of-plane anisotropy, typi- 
cal for Co/Pd multilayers. Sharp antisymmetric peaks 
are clearly observed at about 0.44 T when magnetiza- 
tion reverses its polarity. The antisymmetric peaks are 
superimposed with a slightly asymmetric magnetoresis- 
tance curve. It is important to note that the polarity of 
the odd peaks (positive in the negative field and negative 
in the positive field) is reversed if measurement of the 
longitudinal voltage is done along the opposite edge of 
the film. A similar effect was found by us in FeTb films 
|14j and was previously reported in Co/Pt multilayers 
[S] and (Ga,Mn)As epilayers [7j with perpendicular mag- 
netic anisotropy. Following Cheng et al [6l the odd in field 
longitudinal voltage signal can appear when a domain 
wall separating two domains with up and down magne- 
tization is located between the voltage probes. Electric 
fields generated by the EHE depend on the local magne- 
tization and have opposite polarities on both sides of the 
domain wall. These electric fields normal to the current 
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FIG. 2: (a) Longitudinal voltage Vj (o) measured in the 
Co/Pd multilayer sample at 4.2 K as a function of applied 
field normal to the film plane. Solid line ( ) is a fit according 
to Eq. [s] Dashed lines are a guide for the eye that emphasizes 
the high field asymmetry. Arrows indicate the direction of the 
field sweep, (b) Hall voltages Vti (o) and Vt2 (•) measured 
simultaneously at two locations along the sample. 



can produce a circulating current loop around the domain 
wall and the respective additional voltage contribution 
along the sample. The model of a single domain wall as- 
sumes that magnetization is opposite at locations of the 
two longitudinal voltage contacts when the anomalous 
voltage peaks appear. This assumption can be tested ex- 
perimentally by measuring the Hall voltage at two cross- 
sections along the sample. Fig. |2Jd shows Vt measured 
between two pairs of contacts transversal to the current 
direction, at two locations along the sample (o and •), 
while the longitudinal voltage Vi, shown in Fig. [2^, is 
measured simultaneously between a pair of longitudinal 
contacts. The magnetization reverses almost simultane- 
ously at both locations (the difference in coercive fields 
is about 20 Oe, whereas the reversal width defined as the 
field span over which Vt varies between 10% and 90% is 
approximately 700 Oe (Fig. [2]d)). The antisymmetric 
longitudinal voltage peaks (Fig. [2^) appear with the re- 
versal of magnetization. The width of the peaks is equal 
to the width of the magnetization (Hall voltage) reversal. 
This observation does not agree with the "single domain 
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FIG. 3: Effective circuit representation of the sample. 



wall" picture that predicts opposite Hall voltage polari- 
ties at two cross-sections when the anomalous peaks ap- 
pear. Two more experimental results are important for 
future discussion: (i) Vt signals at two cross-sections are 
similar in shape but differ in magnitude in the magneti- 
cally saturated state at high fields by approximately 8%, 
and (ii) macroscopic magnetization is not uniform: there 
is a small but finite difference in coercive fields along the 
sample. 

Although the single wall model is not in agreement 
with the experimental data, one can assume that the 
transverse voltage is not uniform along the sample. We, 
therefore, model the sample as a simple circuit shown in 
Fig. [3] Vab{H) and Vcd{H) represent the transverse 
voltage generated by Hall effects at two cross-sections 
AB and CD, while two equal resistors R are positioned 
between A and C, and B and D. Field dependence of the 
resistors R{H) is the usual symmetric in magnetic in- 
duction magnetoresistance. Following Kirchhoff 's circuit 
laws, the longitudinal voltages at two edges of the sample 



(2) 



Voltage measured along the sample would differ from 
the ordinary Ohmic if (Vab — Vcd) 7^ 0, i.e. Vab 7^ Vcd- 
The field symmetry of {Vab — Vcd) is even in the case 
of e.g. non-uniform planar Hall effect contribution |15j 
or odd when the ordinary and / or extraordinary Hall 
effects are present. In this case {Vab — Vcd) is given by: 



Vab - Vcd = fJ-o-I 



( Rehe.abMab Rehe,cdMcd 



tAB 



CD 



with tAB and tcD the local thickness, Mab f^nd 
McD as the local magnetizations and Rehe,ab and 
Rem E, CD are the EHE coefficients at cross-sections AB 
and CD respectively. 

Several mechanisms can cause a non-uniform trans- 
verse voltage along the sample. The simplest is a grad- 



ual variation of thickness t ([ij due to either unintended 
imperfection of fabrication or when wedge samples are 
studied. This argument is applicable to any material 
including non-magnetic metals |16j and semiconductors 
|17j . In magnetic materials there are additional mech- 
anisms that can affect the EHE coefficient Rehe- In 
thin ferromagnetic films Rehe depends on the thickness 
and diverges in the thin film limit due to an enhanced 
surface scattering [2 [H]. In granular ferromagnetic or 
superparamagnetic films Rehe depends on size, density 
and shape of magnetic clusters that might not be uniform 
along the sample, due to deposition and annealing pro- 
cedures [4 . If magnetization is uniform along the sample 
{Mab = Mcd), Eq. Ogives: 

- Vcz. = f 1 - f^^^^^^^^ V V4B (4) 

V nEHE.AB ■ tcD J 

The longitudinal voltage Vbd can then be presented 



I R 

Vbd = ^ a ■ Vab 



(5) 



where a is a coefficient that depends on thickness and 
Rehe variation along the sample. The first term on the 
right hand side of Eq. |5] is even with respect to the field 
while the second is odd and proportional to the trans- 
verse voltage. Vab is a monotonic function of field (see 
Fig. [2]d), therefore Eq. [s] can explain the high field asym- 
metry of the longitudinal voltage in Fig. [2^, but not the 
antisymmetric peaks at the magnetization reversal. We 
then assume that magnetization is not uniform and re- 
verses gradually along the sample with raise of the ap- 
plied field. The local magnetization values Mab{H) and 
Mcd{H) are connected by: 

McD {H) = Mab {H) - AH ■ ^^f^"^ (6) 

where Ai/ is the increase of applied field needed to 
propagate the magnetization reversal from the cross- 
section AB to CD. Then, to the first order of AiJ: 



Vac = \ (lR{H) + ^""•''^"•f ^^^'^^ • ^MjAMV 

-1/ if TT}lTT\ ^H-I-flfyREHE.AB dMAB(H 

Vbd = 2 U^l^) 1 mr- 



(7) 



where we assume that Iab — tcD ^ t and Rehe,ab = 
Rehe.cd- The second term in Eq. [Tjis odd with respect 
to the field { H is odd) and can be significant in mate- 
rials with a large EHE coefficient and sharp reversal of 
magnetization, as in thin ferromagnetic films with the 
out-of-plane anisotropy. The shape of AH ■ ^ ^'^'ab(h) ^ 
has a strong peak at magnetization reversal; therefore, 
this term can account for the antisymmetric peaks, as 
in Fig. [2^, or for a significant difference in the maxi- 
mal resistance in Fig. [IJl. If both the gradual reversal 
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of magnetization and the variation of the saturated high 
field Hall voltage along the sample are considered, the 
combination of Eqs. [T] [5] and |7] gives: 



BD 



\ {lR{H) - AH 



dH 



a-VAB{H) (8) 



We applied Eq. j8]to fit the experimental data both for 
the Ni film (Fig. ^) and Co/Pd multilayer (Fig. |2ji) by 
using the measured transverse voltage Vab{H) and two 
fitting parameters and H. The asymmetric magnetore- 
sistance of Ni (solid line in Fig. [iji) was calculated with 
AH = 2.20e and a = 1.4 • 10"^. The fit for the Co/Pd 
multilayer, shown in Fig. [2^ by a solid line, was calcu- 
lated with a = 2- 10^2 ^nd AH = 240e. This value of 
AH is in good agreement with the measured 22 Oe dif- 
ference in coercive fields between cross-sections AB and 
CD. It should be noted that only minor inhomogeneity 
(a) and non-uniformity of magnetization reversal (AH) 
along the sample are sufficient to generate large anoma- 
lous signals. A possible cause for variation of the coercive 
field along the sample is the thickness variation. Magne- 
tization reversal was reported [6l [19] to propagate along 
wedge shaped samples with thickness variation of a few 
percents only. Other possible causes are variation of sur- 
face roughness and adhesion to the substrate which are 
suspected |20j of inducing a transition from nucleation 
dominated reversal to domain-wall-motion reversal. 

Although the model presented above is in a good agree- 
ment with the experimental data, one can wonder if a 
simple circuit (Fig. [s]), which has only two current chan- 
nels, provide a reliable description of a macroscopic sam- 
ple. In the following we present a more rigorous deriva- 
tion of the electric potential along an infinitely long sam- 
ple with variable thickness and Hall resistivity, and show 
that in the proper limit the result is identical to Eq. [8] 
In order to reduce the problem to two dimensions, we 
follow Ref. 16 and define the following two-dimensional 
fields: 



<j{x,y)>=- j{x,y,z')dz' (9) 

^ Jo 



< E{x,y) >E 



t{x,y) 



E{x,y,z')dz' (10) 



where E is electric field, j - the current density, t - the 
average sample thickness and t(x, y) is the actual sample 
thickness at each point. The two dimensional current 
distribution is determined by: 



<E{x,y)>^^-p^-<]{x,y)> (11) 



t{x,y) 



V- <jix,y) >=0 



(12) 



and to a good approximation by: 



Vx < E{x,y) 



(13) 



where ^ is the spatially dependent resistivity tensor. 
Boundary conditions are set to prevent current fiow nor- 
mal to the sample edges. Exponential variation of thick- 
ness along the sample with constant Hall resistivity was 
analysed by Bruls et al [16]. The equations were found to 
be identical to those describing a sample with an expo- 
nential variation of charge carrier density \17\ . Adapta- 
tion of the latter case gives the field dependent potential 
along the sample as: 



'Pix,y,H) 



^p{H) ■ I ■ f3{H) ■ exp (- 



t ■ 



1 



exp 



(14) 



where p{H) is resistivity, a - the length scale over which 
the sample thickness changes by a factor of e, w - the sam- 
ple width, and f3{H) - the ratio of Hall and longitudinal 
resistivities. In the limit of a » x and a » [3{H)w, 
Eq. [14] can be reduced to: 



^{x,y,H) = - 

I p(H)-I-x 



P{H)-I 
t-w 

■_ I p{Hyi-l3{H)-x-{w/2-y) 



x + P{H)-{w/2-y)] 



(15) 



where a — with L being the distance between the 
longitudinal voltage probes and At the change of sam- 
ple thickness between location of the longitudinal probes 
(at x = L/2). The first term in Eq. 15 consists of the 
standard longitudinal and transverse voltages of a homo- 
geneous sample. The second and third terms are correc- 
tions to the potential due to the thickness variation. The 
second term does not contribute to longitudinal voltage 
since it is symmetric in x. The third term is propor- 
tional to the Hall voltage and changes sign depending on 
the location of the probes {at y — Q or y — w). 

Linear variation of Hall resistivity due to change of 
charge carrier density along the sample was analysed by 
Ilan et al [5T] in 2-D electron gas. In the case of mag- 
netic materials we ascribe the gradient of Hall resistivity 
to linear variation of both magnetization and RsHEit) 
along the sample, so that: 



dp^y{x,H) _ 
dx - Mo 



A/(0, H) 
-po ■ Rehe{0) 



ORe 



dt 

dM(x,H) 
dx 



At 
L 



(16) 



where M{0,H) and Rehe{0) are the values of magne- 
tization and Re HE at a; = 0. For p^y varying along the 
X coordinate only, and » w, the potential along the 
sample is given by: 



(p{x,y,H) 



t-w 



x + f3{H) 

p{H)-I-x-(w/2^y) 



{w/2 ~ y)] 



(17) 
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where 1^{H) = p{H)-(^ 



dx 



211. The first term 



in Eq. [17] corresponds to the potential distribution in 
a homogeneous sample, whereas the second term is the 
correction due to a spatial variation of the Hall resistivity. 
Since the correction terms in Eqs. 15 and [17] are small 



and of different origins, they are additive (higher order 
corrections are neglected). The longitudinal voltage can 
be calculated from Eqs. [15] [T7|as: 



Vi{H)^ip{x^-L/2,H) 



At 



I /■MO-M(0,g)-flEHE(0) 

2-t \Rehe{0) 

L-J-Mo--R.Effg(0) dM{x,H 
2-t ' dx 



L/2,H) - 

dREHEit) _ Ai' 
dt t 



±- 



where the ± sign stands for y — (+) and y 
The first term in Eq. 



(18) 
w{-). 



18 is simply IR{H), the second 
term is the correction due to the thickness variation and 
is proportional to the transverse voltage, and the last 
term is the correction due to a non-uniform magnetiza- 
tion along the sample. Finally, by assuming a constant 
ratio between the change of an applied field AiJ and the 
propagation of the magnetization reversal over a distance 
L, we calculate: 



Vi{H)^I-R{H)±\(^ 



At 



± 



^ dRE 

Rehe{0) 
AH dVt(0,H) 



dt 



'At) 



dH 



(19) 
is identi- 
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with VtiO,H) = i-^o-REHE{oyAmH) 
cal to Eq. [8| obtained from the circuit modelT 

The analytic calculation was done for an infinitely long 
sample. In order to treat a finite sample, numerical cal- 
culations were carried out. Following Hajjar et al [25], 
the current distribution was calculated by taking the 
finite difference version of Eqs. [12] - Jl3] on a two di- 
mensional rectangular lattice. Boundary conditions were 
added along the length of the sample edges together with 
the current source and drain. The sample dimensions 
were chosen equal to the actual geometry of the Co/Pd 
multilayer film, in which L = w and the total sample 
length is 3w. The effective resistivity tensor that was 
used included both thickness and magnetization gradi- 
ents along the sample: 



t-p{x,H) 
t{x) 



1 



At 
t-L 



P{H) ( 1 + 

-PW - ^ 1 

(20) 

In order to obtain lx{H) at each field value, the nor- 
malized change in magnetization between locations AB 
and CD was estimated as: 



AM{H) Vcd{H) Vab{H) 



Mo 



Vc 



CD.S 



Va 



B,S 



(21) 



where Ms is the saturation magnetization, Vab.s f^nd 
VcD,s are the saturated values of Vab and Vcd respec- 
tively. The normalized magnetization slope along the 
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FIG. 4: Longitudinal voltage Vi (o) measured in the Co/Pd 
multilayer sample at 4.2 K as a function of field normal to the 
film plane in the peak region. Solid line (-) is a fit calculated 
by using Eqs. 



E 



20 21 dashed line (- -) is a fit according to Eq. 



sample is then ^^-^j^- The solid line in Fig. jZ] presents 
the simulation of the field dependent longitudinal voltage 
in the peak region for the Co/Pd multilayer sample with 
a single fitting parameter Ai/t = 0.05. The dashed line 
was calculated by Eq. J8] and is shown here for compari- 
son. Numerical results agree nicely with the experimental 
data (o). A snapshot of the simulated sample potential 
during a gradual magnetization reversal (magnetization 
is zero at x = at applied field of 0.44 T) is shown in 
Fig. 5 Large Ohmic component was subtracted 

for clarity. It is clearly seen that the potential gradi- 
ent along the sample has opposite polarities at two edges 
of the sample. It is important to note that due to the 
sharpness of magnetization reversal in films with an out- 
of-plane anisotropy a minor delay in coercive field (20 Oe 
as compared with 700 Oe of the reversal width) results in 
a relative difference of magnetization of up to about 20% 
between cross- sections AB and CD, which respectively 
leads to distinctive voltage peaks in magnetoresistance. 

To summarize, we studied the asymmetric field de- 
pendence of magnetoresistance in magnetic films. We 
argue that minor variation of thickness. Hall coefficient 
and nonuniform magnetization reversal along the sam- 
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pie can explain the anomalous phenomena. We show 
that a non-uniform variation of the Hall voltage along 
the sample generates an additional odd in field longitu- 
dinal voltage signal proportional to the field derivative 
of the transverse voltage. This additional signal can be 
significant when the Hall voltage varies sharply with the 
applied field, like in the case of magnetization reversal 
in films with perpendicular magnetic anisotropy, stud- 
ied here, at superconducting transitions or in materials 
demonstrating the quantum Hall effect. The fingerprint 
of the mechanism is the reversal of the asymmetry when 
the longitudinal voltage is measured along the opposite 
edge of the sample. 



FIG. 5: Numerical calculation of electric potential generated 
by a nonuniform magnetization reversal. The standard I ■ R 
contribution is subtracted for clarity. A,B,C and D corre- 
spond to locations of the voltage probes. 
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